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1. Introduction. 



New insights into the behaviour of low energy supersymmetric gauge theories 
have led to an explosion of activity in the subject [1]. One aspect which has at- 
tracted attention is the low-energy properties of the Coulomb phase of N=2 gauge 
theories with Nf matter multiplets in the fundamental representation, with and 
without bare masses. Since the moduli space of these gauge theories coincides 
with that of a particular hyperelliptic curve [1]— [15], [23, 24], one attempts to ex- 
tract the strong-coupling physics from the curve describing the theory in question. 
The desired information can be found from the prepotential jF(a) characterising 
the low-energy effective Lagrangian, and from the Seiberg-Witten (SW) period 
integrals 

— * 

n = 

These are related to J 7 (a) by 

In the weak-coupling region one can write (schematically) 

F{a) = T i_i oop (a) + •^instantons( a )- ( L3 ) 

A number of different methods have emerged for evaluating equation (1.3) for a 
variety of groups, either using the properties of hyperelliptic curves, or from the 
microscopic Lagrangian itself [30]. However, the available results for the strong- 
coupling regions are still rather sparse, and as yet no complete picture has emerged, 
except for gauge groups of rank 1 or 2. 

One strategy for obtaining strong-coupling information is to derive and solve a 
set of Picard-Fuchs (PF) equations for the SW period integrals. The PF equations 
have been derived in a number of special cases for Nf = and massless multiplets 
for N f ^ [8, 19, 20, 21, 29]. 
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The PF equations have also been formulated [20, 21, 29] as an explicit, analytic 
set of PF equations valid for arbitrary classical gauge groups, and an extensive 
set of values for Nf massless multiplets in the fundamental representation. A 
systematic method for finding PF equations for cases with bare mass m — 0, which 
is particularly convenient for symbolic computer computations, was given by Isidro 
et al [21]. Alternative techniques, not involving PF equations, are explored in [16, 
17]. 

The situation for m ^ is less fully developed [1, 18, 22, 27, 28]. Part of the 
problem is that the preferred SW differential Xsw, is an abelian differential of the 
third kind if m ^ 0, while of the second-kind if m — 0. 

The SW differential has the property that 




where 7$ and 7^ are closed 1-cycles with canonical intersection in a homology 
basis. Since Xsw is of the third- kind when m 7^ 0, i.e., it has poles with non- 
vanishing residues, some care is required to formulate the PF equations. Explicit 
PF equations for SU(2) and SU(3) have been presented for m/fl [22, 28]. Even 
for these simplest examples several pages are required to write out the actual PF 
equations. Thus, it does not seem practical to present explicit PF equations for 
arbitrary gauge groups and Nf consistent with asymptotic freedom when m 7^ 
0. Rather, it seems more convenient to present a comprehensive set of explicit 
algorithms, from which one may obtain the PF equations with m^0. 

In this paper we generalise the results of [21], and present a new method to 
obtain the PF equations for arbitrary classical gauge groups, with massive matter 
hypermultiplets in the fundamental representation. We have obtained the explicit 
PF equations for a number of gauge groups with massive multiplets. The corre- 
sponding PF equations turn out to be extremely lengthy. Therefore, there is no 
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virtue in presenting these results. However, explicit results using our methods co- 
incide with PF equations previously presented in the literature for m 7^ 0. Thus 
our paper concentrates on a description of our method. In section 2 we formulate 
the problem, while in section 3 a necessary set of recursion relations are exhibited. 
Sections 4 and 5 describe how the PF equations are derived from this information. 
Final comments, as well as a summary and conclusions, appear in sections 6 and 
7, respectively. 



2. Formulation of the problem. 

The strategy used in [21] to derive the PF equations of effective N = 2 super- 
symmetric Yang-Mills theories in 4 dimensions can be modified to include massive 
matter hypermultiplets. To do so, let us first recall the necessary elements from 
[21]. 

Let us consider an effective N = 2 supersymmetric Yang-Mills theory charac- 
terised by a certain gauge group G with rank r, a number Nf of massive matter 
hypermultiplets with bare masses mj, where 1 < j < Nf, and a number of moduli 
Ui, where 1 < % < r. The matter hypermultiplets will be taken in the fundamental 
representation, and Nf will be restricted to those values for which the theory is 
asymptotically free, but for the moment G will remain unspecified. Consider the 
complex algebraic curves 

y 2 = P(x;ui; mj ;A), (2.1) 

where P is a polynomial in x of degree 2g + 2, and A is the dynamically generated 
quantum scale. When all roots of P are pairwise different, i.e., away from the 
zero locus of the discriminant of the curve, equation (2.1) defines a family of non- 
singular hyperelliptic Riemann surfaces E 5 of genus g [31]. Under an appropriate 
choice of the polynomial P, the moduli space of quantum vacua of the theory under 
consideration is coincident with that of the curves defined by equation (2.1). On 
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Ylg there are g holomorphic 1-forms 



• n t 

x>—, j = 0,l,...,g-l. (2.2) 

y 

The following g 1-forms are meromorphic on T, g and have pole singularities at 
infinity of order greater than 1 * 

• H t 

x J — , j = g + l,g + 2,...,2g. (2.3) 

y 

Furthermore, the 1-form 

X 9— (2.4) 

y 

is meromorphic on S fl , with a simple pole at infinity. Altogether, the abelian differ- 
entials x J dx/y in equations (2.2) and (2.3) will be denoted collectively by cuj, where 
j — 0,1, ... , 2g, j ^ g. We define the basic range R to be R — {0, 1, . . . , g, . . . 2g}, 
where a check over g means the value g is to be omitted. 

Following [21] we give two definitions. Let us call W = y 2 = P(x;uf, my, A) in 
equation (2.1). Moreover, given any differential x n dx/y, with n > an integer, let 
us define its generalised fi-period (uf, my, A; 7) along a fixed 1 -cycle 7 G Hi(E g ) 
as the line integral [32] 

^(ui-mj-A;^ := (-1)^(^ + 1) f^dx. (2.5) 

7 

In equation (2.5), T stands for Euler's gamma function, while 7 G H\(llg) is any 
closed 1-cycle on the surface. For the sake of simplicity, we will drop m, mj, A and 
7 from the notation for the periods . As explained in [21], we will work with 
an arbitrary value of \i, which will only be set to —1/2 at the very end. 



* The differential x 9+1 dx/y is of the second kind, i.e., it has vanishing residues, for the 
particular curves of Seiberg-Witten theories. The differentials x n dx/y have non- vanishing 
residues at infinity when n > g + 1. This corrects a mistake in the terminology of [21]. 



■5 



In effective N = 2 supersymmetric gauge theories, there exists a preferred 
differential, called the Seiberg-Witten (SW) differential, \sw, with the following 
property [1]: the electric and magnetic masses a« and af entering the BPS mass 
formula are given by the periods of Xsw along some specified closed cycles 7$, G 
Hi(Eg), as in equation (1.4). In these theories, the polynomial P(x; uf, mj; A) is of 
the special form P = p 2 (x) — G(x), for certain p(x) and G(x) (see below). Then, 
the SW differential takes on the following expression: 



The SW differential further enjoys the property that its modular derivatives 
dXsw I 'dui are (linear combinations of the) holomorphic 1-forms [1]. This ensures 
positivity of the Kahler metric on moduli space. 



As seen in the massless case treated in [21], the periods VLn defined in equation 
(2.5) satisfy a set of recursion relations in both indices n and \x that can be used to 
derive the PF equations. Similar conclusions continue to hold for the massive case 
as well. However, the polynomial P{x\ u\\ m,j] A) defining the curve depends on the 
particular gauge group G in such a way that a general expression valid for all G, 
such as the one used in [21], cannot be given. Instead, the different gauge groups 
have to be treated separately. As the derivation of the recursion relations follows 
the same pattern used in [21] for the massless case, it will not be reproduced here. 
We will simply list the final results below. 




(2.6) 



3. The recursion relations. 



a ) G = SU(N C ), N f < N c . 



The curve is given by [23] 



W =p 2 (x)-A 



2N C -N f 



Nf 

Y[(x + mj), 



(3.1) 
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where 



P(x) = ^2uiX l , U Nc = 1, WjV c _i = 0. 



(3.2) 



j=0 



Defining the symmetric polynomials in the masses SN f -j(m) through the expansion 

(3.3) 



Y[(x + mj) = J2 S N f -j(m)x J , 



one finds that the following recursion relations hold: 



iV c -ljV c -l iVc-1 

E E»« !i ^-2E^ ( " +1) 



w 

n 



(3.4) 



iVc+n+j 



and 



i=0 i=0 



-f 



j=0 



E 

J=0 



A£--"'£(aw«-i)s w ,_,i#/> 



n + 1 - 2JV C (1 + //) 

7V C -17V C -1 7V C -1 

+ E E 0' + « - 2N c )« i « J n^ + 2 5> - tf c )«;njS£^ 



(3.5) 



j=0 Z=0 j=0 

When n + 1 — 2 iV c (2 + //) ^ 0, one can combine equations (3.4) and (3.5) to obtain, 
after shifting n + 2N C — > n, 

1 



q(m+i) 



n + 1 - 2iV c (2 + /i) 



x 



3=0 



N c -1 



+ 2 £ ((i + ^)(iV c + J V(n-2Ar c + l)) M X- + A ) +J 



(3.6) 



3=0 
N C -1N C -1 

+ E E (0" + l)(l + l*)-(n- 2N C + l)) Ui u,n<£^_ 



2N C 



Modular derivatives of periods are given by 



E 

3=0 



^ = 25>«£& (3.7) 



b ) G = SU(N C ), 2 < N c < N f < 2N C . 
The curve is given by [23] 



W 



Nf—N c N f 



3=0 3=1 



with as in equation (3.2). We observe that this curve can be obtained from 
the one given in equations (3.1) and (3.2) by a shift of some of the moduli uf 



u t ^ Ui + ^Af f c Nf S Nf - Nc _i, 0<i<N f -N c . (3.9) 



In the classical limit, this redefinition does not affect the moduli. Applying this 
shift to equations (3.4) through (3.7) we can straightforwardly derive the recursion 
relations corresponding to this case from those of the previous case. Alternatively, 
one could carry out a step-by-step derivation like the one needed for Nf < N c . 
Either way, the recursion relations turn out to be given by 



Nf N C -1N C -1 
Q (H+1) _ .2N C -N f y- „ V V 7/ -7/,0 (/i+1) 

* l n+2N c - A N f 2^ bN f-r L n+j ~ 2^ U 3 U ^ l n+j+l 

3=0 j=Q 1=0 

N c -1 N c Nf—N c 

- 2 e - \<r N ' E E ^<l + „- i+1 (3.10) 

j=0 1=0 i=0 



Nf—N c Nf—N c 

,(M) 

16 

3=0 1=0 
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and 



n + 1- 2N C (1 + fi) 



N f 
3=0 



(3.11) 



N c -1 N C -1N C -1 

+ 2j2u- n c )u^i +j + E E 0' + ' - 2N c)«i«ito!t£l 

3=0 j=0 1=0 

Nf—N c Nf—N c 

+ s A " E E ( 2JV / - 4Ar = - j - wm'£$,,- w .- i -i 

3=0 1=0 

N c Nf—N c 

z=o i=o 

When n + 1 — 2iV c (2 + //) ^ 0, one can combine equations (3.10) and (3.11) to 
obtain, after shifting n + 2N C — > n, 



TV 



n + 1 -2iV c (2 + /x) 



< e ""' E (» " 2N < + 1 - ^ + »))S Nf -3^ Nc+j 

3=0 



pAN c -2N s Nf—N c Nf—N c 

— - E E (n-2iV c +l- (1 + ^(2^ -2JV c -j-0)^n^ / _ 4JVe _ J ._ J 



16 

z=o j=o 



N c Nf-Nc 



iV c iV/-iV c 

E E 

1=0 j=o 

_l 1 _ ft _l _l ;^„, .„,.o(m- 



AT C -1 JV C -1 

]T ]T (n - 27V C + 1 - (1 + M + Z)) Uj .„ z n^S} 
j=0 z=o 
AT C -1 

2 £ (n - 2iV c + 1 - (1 + + N c )) Uj rt n ^l +j 

j=0 



(3.12) 



Modular derivatives of periods are given by 

Nc - 

~ A 2-^ U ^ L n+i+j + n IV N f 2^ a ^ l N f -N c +n-j+i- y 6 ' 16 ) 



c ) SO(N c ),N f <N c -2. 



The curve for Nf < N c — r — 2 is given by [24] 



N f 

W = p 2 (x) - k 2 ^- Nf ~ 2) x d JI^ 2 - ™ 2 ), (3.14) 



with the rank r and the power d being respectively given by r = N c /2 and d — 4, 
if N c is even, and r = (N c — l)/2 and d = 2, if iV c is odd. The polynomial is 
given by 

= ^MjX- 7 , U 2r = 1, Modd = 0. (3.15) 

i=o 

Let us expand the mass term of equation (3.14) in terms of the symmetric polyno- 
mials in the squared masses TN f ~j{m 2 ) as follows: 

N f N f 

l[(x 2 - m 2 ) = J2(-l) Nf - j TNf-j(m 2 )x 2 i. (3.16) 

j=l j=0 

The above expansion could just as well be expressed as a double summation with 
coefficients given by the symmetric polynomials in the masses Sj\f f ^j(m), as done 
for SU(N C ) in equation (3.3). However, we will find it more convenient to use the 
expansion (3.16), as it manifestly preserves the even parity of the SO(N c ) curve 
under x — > —x. Reasoning as in [21], one finds that the following recursion relations 
hold: 



N 



n+1 — 4r(l + ji) 



n+d+2i 

i=0 

2r-2 2r-2 2r-2 

(P+ 1 ) , o V (a o^„. .o(a*+ 1 ) 



+ E E (* + j ~ + 2 E V - 2r)^^ 2 . 

i=0 j=0 j=0 



(3.17) 
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and 



(M+l) 



1 



n + 1 - 4r(2 + /i) 

at, 



■x 



A 2(N c -N f -2) J2(-l) N ^(n - 4r + 1 - (1 + rf(2i + d))T Nf ^ n+d _ 4r+2i 



i=0 
2r-2 

+ 2 £ ((1 + ^(j + 2r) - (n - 4r + l))^- + 2r +j 

3=0 
2r-2 2r-2 

i=0 j=0 



Modular derivatives of periods are given by 



j'=o 



(3.18) 



(3.19) 



It is known from [24] that the curve for N c — r — 2 < Nf < N c — 2 can be 
obtained from that for Nf < N c — r — 2 by a shift of the moduli similar to the 
one performed in equation (3.9). The necessary shift now affects the even moduli 
only, and it involves the symmetric polynomials in the squared masses Tj(m 2 ) 
rather than the Sj(m). Application of this shift to the recursion relations given 
in equations (3.17), (3.18) and (3.19) will produce the recursions corresponding to 
N c - r - 2 < N f < N c - 2. 

d ) Sp(N c ), < N f < N c + 2. 

The curve is given by [5, 25]* 



x 2 W 



p 2 (x] 



-A^ +2 -^n(- 2 + 

3=1 



(3.20) 



★ The sign of the m 2 term in ref. [25] is opposite to that of eqn. (3.20). If so, that would mean 
the double scaling limit would only apply to an even Nf. The replacement TM f -j{m 2 ) — ► 
(_ l} N f~iT N j(m 2 ) in cqns. (3.22), (3.23) and (3.24) below will accommodate the sign 
given in ref. [25]. 
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where N c = 2r. The polynomial p(x) is given by 

2r+2 N f 

p( x ) = e u j x ° + A N^ 2 ~ Nf n u2r + 2 = i? M ° dd = °- ( 3 - 21 ) 

J'=2 J=l 

The mass term of equation (3.20) can be expanded in terms of the symmetric 
polynomials Tjsr f -j{m 2 ) defined by 



N f 



(3.22) 



j=l j=Q 

Similar steps to the ones taken above then lead to the following recursion relations: 



1 



n+l-(4r + 2)(l + /i) 



Nf 



2r 2r 2r 

+ 2 (i ~ 2(r + l)H4+ + 2^ + EE( i+ ^ 4 ( r + l))^«i*W;-2 

j=2 j=2 



Oi+i) 

n+2j-2 



i=2 

iV/ 2r+2 

i=l j=2 



(3.23) 



and 



q (m+i) 



2r 2r 



n + 1- (4r + 2)(2 + /x) 



x 



[ £ £ ((i + , - 2)(1 + - (n - 4r - l))^^<^-4,-4 



i=2 j=2 
2r 



+ 2 J3 ((1 + /i)(2r + j) - ( n - 4r - 1))^ 



Ou+l) 
n-2r+i-2 



i=2 



JV/ 2r+2 



(3.24) 



[n 



+ 2A^ +2 -^n^E(( i +^'- 2 ) 

i=l j=2 

+ A?^ 2 -^ E «» - 4r - 1) + (1 + „)(2 - 2j))T NM ^l 



4r-l))«inJS-4r-4 



-4r-4 
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Modular derivatives of periods are given by 



ant/i) 2r+2 N f 

^ = 2 £ u^_ 2 + 2A™-»> J] (3.25) 

j=2 j=l 



4. Derivation of the Picard-Fuchs equations. 

Following the steps of [21], one can use the recursion relations of the previous 
section to derive a coupled system of first-order, partial differential equations (with 
respect to the moduli) satisfied by the periods. We first set /i — —1/2 in all what 
follows. Then we need to identify the appropriate subspace of periods that one 
must restrict to, in order to properly solve the above recursions. We recall that a 
key element is the behaviour of the curves under the operation of parity x — > —x. A 
glance at the equations of the previous section immediately reveals that conclusions 
completely analogous to those of [21] continue to hold for the massive case. Let us 
briefly recall them. 

• For the SO(N c ) gauge groups, one must restrict to the even subspace of the 
basic range R, i.e., to even values of the subindex n. This follows from two facts. 
One is that all the odd Casimirs of the gauge group SO(N c ) vanish, so the recursion 
relations have a step of 2 units. Furthermore, the solution to those recursions can 
be expressed in terms of a set of initial data with an even value for the subindex. 
As a function of the rank r, the genus g of the SO(N c ) curve is g — 2r — 1. This 
being odd, the value of the subindex n at which the recursion (3.17) blows up is 
skipped. Similarly, the zero of the denominator of equation (3.18) is avoided when 
n is even. 

• For the Sp(N c ) gauge groups, one must restrict to the odd subspace of the 
basic range R, i.e., to odd values of the subindex n. Again, this follows from the 
same facts as above. All the odd Casimirs of the gauge group Sp(N c ) vanish, so the 
recursion relations have a step of 2 units, but the solution to those recursions can 
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be expressed in terms of a set of initial data with an odd value for the subindex. 
As a function of the rank r, the genus g of the Sp(N c ) curve is g — 2r. This being 
even, the value of the subindex n at which the recursion (3.23) blows up is skipped. 
Similarly, the zero of the denominator of equation (3.24) is avoided when n is odd. 
The factor of x 2 present in the left-hand side of equation (3.20) is responsible for 
this odd parity, as opposed to the even parity of the SO(N c ) recursions. 

• For the SU (N c ) gauge groups, the curve has no well defined parity under 
x — > —x. This is a consequence of the fact that SU(N C ) has both even and odd 
Casimirs which, in turn, causes the recursion relations have a step of 1 unit. To 
solve the recursions, one must first work on the enlarged subspace of periods given 
by RU{g}, i.e., fit ±1/2) , with n E RU{g}. This is in order to avoid the divergence 
of the recursion relations (3.5) and (3.11) that occurs when the subindex n takes on 
the value n — g, where the genus g now equals the rank r, g = r. Next, one applies 
a linear relation satisfied by the fi^ +1 ^ 2 \ where n G R U {<?}. For a derivation of 
this linear relation that also holds in the massive case treated here, see [21]. 

Once the correct subspace of periods has been identified, the recursion re- 
lations can be solved as explained in [21]. Let us arrange the periods fil ±1//2 ' ) , 
where n spans the appropriate subspace of R, as column vectors: f2( ±1//2 ) = 
(fij ±1/2) , . . . , ftJ ±1/2) , ^+ 1 / 2) , . . . , ^ 1/2) Y- We have called the dimension of the 
subspace under consideration 2r, as it turns out to equal two times the rank r 
of the gauge group. We have also arranged the entries of fi^ 1 / 2 ) in such a way 
that the first r of them are holomorphic, while the last r of them are meromorphic. 
From the recursion relations, the vectors of periods f^ -1 / 2 ) and f2( +1 / 2 ) are linearly 
related through a matrix M 

fi(- 1 /2) =M . fi (+l/2) ? 



where M is (2r x 2r)-dimensional. Its entries are certain polynomial functions in 
the moduli U{, the bare masses mj and the quantum scale A, explicitly computable 
using the recursion relations given in the previous section. 
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Similarly, one can use the expressions for the modular derivatives (equations 
(3.7), (3.13), (3.19) and (3.25)), together with the recursions, in order to write a 
system of equations which, in matrix form, reads 

A n (- 1 /2) = £,(„.). n (+ 1 /2). (4.2) 

OUi 

The matrix D{ui) is (2r x 2r)-dimensional. Again, its entries are some polynomial 
functions in the moduli Ui, the bare masses mj and the quantum scale A, explicitly 
computable from the recursion relations. Assume for the moment that the matrix 
M of equation (4.1) is invertible. Combining the latter with equation (4.2) one 
ends up with 

A n (-i/2) = v . . n (-i/2) (4 . 3) 

dm 

where we have defined the matrix Ui as 

Ui = D( Ui ) ■ M-\ (4.4) 

Equation (4.3) is a coupled system of first-order, partial differential equations satis- 
fied by the periods fi( -1 / 2 ): the first-order PF equations. It expresses the modular 
derivatives of the basic periods fi( -1 / 2 ) (of the subspace under consideration) as 
certain linear combinations of the same periods Q( -1 / 2 ). The coefficients entering 
those combinations (i.e., the entries of the Ui matrices) are some rational functions 
of the moduli Ui, the masses rrij and the quantum scale A, explicitly computable 
from the recursion relations. 

To close this section, let us comment on the invertibility of the matrix M of 
equation (4.1). Explicit evaluation in a wide class of examples shows det M to 
be a product of the factors of the discriminant A(ui,mj, A) of the corresponding 
curve, possibly with different multiplicities, and up to an overall non-zero constant. 
Therefore, the M matrix encodes the singularity structure of the curve, and it is 
invertible except at the singularities of moduli space, i.e., except on the zero locus of 
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the discriminant A(u{, rrij, A) of the curve. This conclusion holds with some caveat 
when the gauge group is SU(N C ), since det M may then pick up some additional 
zeroes. This possible new zero locus of det M occurs for the same reasons already 
explained in [21] for the massless case. The inclusion of non-zero bare masses rrij 
does not alter the arguments given in [21] as far as det M is concerned, and they 
continue to hold in the massive case as well. 



In principle, integration of the system (4.3) yields the periods as functions of 
the moduli U{. The particular 1-cycle 7 e ifi(E fl ) being integrated over appears in 
the specific choice of boundary conditions that one makes. In practice, however, 
the fact that the system (4.3) is coupled makes it very difficult to solve. A possible 
strategy is to concentrate on one particular subset of periods and try to obtain a 
reduced system of equations satisfied by it, at the cost of increasing their order. 
In [21] we have made use of the fact that one can perform a change of basis 
that included Slsw-i the period of \sw-, as a basic vector. The decoupling of the 
resulting equations then followed from the property that the modular derivatives 
of the SW differential \$w are the holomorphic differentials of the appropriate 
subspace within which the recursions are being solved — we call this the potential 
property of Xsw- However, a similar change of basis is inconvenient now, because 
Xsw is °f the third kind. Let us see how this difficulty can be circumvented. 

Consider the U{ matrix in equation (4.3) and block-decompose it as 



where all four blocks B{, C{ and D{ are r x r. Next take the equations for the 
derivatives of the holomorphic periods, dfl n /dui , 1 < n < r, and solve them for 
the meromorphic periods Q n , r < n < 2r, in terms of the holomorphic ones and 



5. Decoupling the Picard-Fuchs equations. 




(5.1) 
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their modular derivatives. That is, consider' 



_d_ 

dm 









/ Q r+ i \ 
















I ^2r / 



(5.2) 



Solving equation (5.2) for the meromorphic periods involves inverting the matrix 
B{. Although we lack a formal proof that B{ is invertible, when the rank r of the 
gauge group is greater than 1, det Bi turns out to have two types of zeroes on 
moduli space. The first zero locus contains a product of the factors of the discrim- 
inant A(ui,m,j,A), possibly with different multiplicities. The second zero locus is 
unrelated to the discriminant. Away from these singularities, B{ is invertible so, 
from equation (5.2), 



V ^2r / 



(5.3) 



Next, substitute the meromorphic periods (5.3) into the last r equations of (4.3), 

A ; =q : +A : i. (5.4) 



/ Q r+ 1 \ 
















+ A 




\ n 2r ) 








V ^2r / 



and rearrange terms in order to obtain 



t 



+ 



\n r J 



dB t 



dui 1 



B7 l + Ai + BiDiBr 1 



_d_ 

dui 



Ihlhll '.1; IhC; I ^H. '.1; -° A> 



dm 1 dm 



(5.5) 



Equation (5.5) is a second-order coupled system, satisfied by the holomorphic pe- 



* For notational simplicity we have dropped the superscript /i = —1/2, with the understanding 

that it has been fixed, 
t No summation over i is implied here. 
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riods on the curve: the second-order PF equations. In order to decouple them, one 
now employs the potential property of the SW differential Xsw, 



_d_ 

dm 



Xsw = ^i, i = l,...,r, 



(5.6) 



and from here one concludes an analogous property for the corresponding periods, 



d 

OUi 



(5.7) 



The passage from equation (5.6) to equation (5.7) is justified, even though in the 
presence of non-zero bare masses mj the SW differential Xsw is of the third kind. 
The reason is that, from [1], the residues of Xsw are known to be some multiples 
of the bare masses rrij, and are therefore independent of the moduli U{. 

Finally, substitution of equation (5.7) into (5.5) produces a decoupled system 
of third-order, partial differential equations for the SW period Slsw, 



+ 



/ diQsw \ 



\ d r tt S W I 



OB. 



t Br 1 + A i + B i D i Br 1 



du- 1 





( di^sw \ 


d 


dui 






\ o r n S w J 



dBi 



BiDiB^Ai - Bid + g^B^A 



dm 



( di^sw \ 



(5.8) 



0, 



\ d r flsw I 

which are the third-order PF equations of the massive N = 2 theory. 
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6. Final Comments. 



We observe that the recursion relations on which our method is founded can 
be derived without taking recourse to the SW differential \$w or its period QsW- 
all that is required is a knowledge of the curve. One first derives a set of recursion 
relations. Next one solves them, to obtain a first-order system (4.3) whose sub- 
matrices Ai, B{, Ci and Di contain all the relevant information. The final step is 
the decoupling procedure. Therefore, whatever limits we may want to take in the 
equations above must be taken at the level of the recursion relations. If the latter 
enjoy the correct limiting properties, so will the PF equations derived from them. 

Two limits are worth taking. One is that in which all bare masses rrij tend to 
zero: rrij — > 0. This is the reduction to the massless case. The other one is the 
integrating out of one massive quark, also called double-scaling limit: sending the 
j-th mass to infinity, rrij — > oo, and the quantum scale to zero, — > 0, while 
keeping rrijl\ p Nf constant and setting it equal to the new quantum scale A p Nf _ v 
(j) and p' being the required powers to which A^ f and A^-i are raised, in the 
presence of Nf and Nf — 1 flavours, respectively). This removes one flavour from 
the problem. One can easily check that all the recursion relations given in section 
3 tend to the corresponding recursions (either massless, or with one massive quark 
less) in the appropriate limits. This is a trivial consequence of the fact that the 
massive curves themselves are so constructed as to reproduce both limits correctly. 
For the massless limit in particular, the recursion relations already derived in [21] 
are also correctly reproduced. 

However, the limit in which all quarks become massless is more intriguing, in 
the following sense: the final equations for flsw are third-order in the massive case, 
while they are only second-order in the massless case. How does this reduction in 
the order of the equations take place? This point has been argued in [22], but let 
us see how it can be recast in our language. 

As already remarked, the system of first-order equations (4.3) not only enjoys 
the correct limiting properties, but its derivation follows the same pattern in both 
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the massive and the massless case, since no use is made of the SW differential Xsw 
or its period Qsw- However, the passage from the first-order equations (4.3) to 
the second-order equations (5.5) (or equations (3.10) of [21]) differs in the massive 
and the massless cases. In both cases one solves for the meromorphic periods in 
terms of the holomorphic periods and their modular derivatives. However, in the 
massless case one can immediately apply the potential property of the SW period 
flsw, while in the massive case one still substitutes all the meromorphic periods 
into the remaining first-order equations. This further step must be taken because, 
the SW differential now being of the third kind, one cannot apply the potential 
property directly. This accounts for the increase in the order of the final equations. 

The limit mj — > in which all masses vanish must therefore be taken in 
the recursion relations, or at most at the level of the first-order equations (4.3), 
i.e., prior to the decoupling of the equations. The reason is that the decoupling 
procedure does not commute with the limit, as it proceeds differently in the two 
cases. 

7. Summary and Conclusions. 

In this paper we have extended a previously described derivation of the Picard- 
Fuchs equations [21], in order to include the case of massive matter hypermultiplets 
in the fundamental representation. It is systematic, well suited for symbolic com- 
puter computations, and holds for any classical gauge group. Our method is based 
on a set of recursion relations satisfied by the period integrals that one can define 
on a hyperelliptic Riemann surface. We explicitly focused on the case of effective 
N = 2 supersymmetric Yang-Mills theories in 4 dimensions, where the relevant 
Riemann surfaces are such that their moduli space coincides with the moduli space 
of quantum vacua of the theory under consideration. From this point of view, the 
Picard-Fuchs equations have proved to be an important tool in probing the struc- 
ture of moduli space and in computing the full prepotential, including instantons. 
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For gauge groups with rank 2 or greater, the Picard-Fuchs equations in the 
presence of massive matter become intractable. The goal of explicitly computing 
the full quantum prepotential of effective N = 2 theories may well have to be 
accomplished using alternative techniques, such as the ones put forward in [16, 17]. 
However, we believe the method presented here may find application elsewhere, as 
our derivation is not limited to these specific areas, and its algebraic nature lends 
itself easily to different uses. 
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APPENDIX 

Computer Calculation of the PF Equations 

Let us return to equation (5.5), where we have observed that the entries of 
the blocks A{, B{, Ci and D{ are certain rational functions of the moduli Ui, the 
bare masses mj, and the quantum scale A. This rational character stems from the 
inversion of the M matrix in equation (4.4), where division by det M is required. 
As explained, det M is also a polynomial in the moduli and the bare masses. 
Transferring detM to the left-hand side of equation (4.3) leaves us with 

det M^(- 1 /2) = i/..Q(-V2) ) (A.l) 
dui 

where Ui is the matrix product D(ui) ■ M _1 , M being the matrix of cofactors of 
M. Given that both D(ui) and M have polynomial entries in the moduli and 
the masses, U{ also has purely polynomial entries, rather than rational functions* 
To the effect of performing symbolic computer calculations, this has an obvious 
computational advantage, at the small expense of having to modify equations (4.3) 
through (5.5). We first block-decompose Ui in equation (A.l) as 

ft -(££)■ (A - 2) 

Next we solve for the meromorphic periods in terms of the holomorphic ones, as 



★ The one possible exception to this polynomial character of Ui may (but need not) occur for 
the SU(N C ) gauge groups, under the circumstances explained in [21]. 
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in equation (5.3), 



/ &<r+l \ 
V ^2r / 



(A.3) 



and substitute into the remaining equations, as done in (5.4). After some rear- 
rangements we find^ 



^-^(detM)- 1 ^ 
dm 



d 

+ (det M^BiDiB' 1 - (det M)" 1 — det M 

OUi 







d 


dui 









(A.4) 



+ 



(det M)- 2 B l b l B^ 1 A i - (det M)~ 2 BiCi 



+ (det My 



-id** 



dm 



R-Mi-(det My 



dm 



/OA 



Obviously, equation (A.4) reduces to (5.5) upon formally setting det M = 1 and 
eliminating the tildes, as this corresponds to leaving det M in the denominator in 
the right-hand side of equation (A.l). 

Below we list a Mathematica programme that computes the matrices Ai, Bi, Ci 
and L>i, as well as det M, for the SO(N c ) gauge theory with Nf massive multiplets, 
when Nf < N c — r — 2. It can be easily modified to apply to the SU(N C ) and 
Sp(N c ) gauge groups. For typographical reasons, the names of some variables in 
the body of the programme have been changed with respect to the paper. Thus, 
P n stands for ^ +1/2) , Q n for fiJr 1/2) , "mu" for /i, "nf for Nf, "nc" for N c , R for 
A(ui,mj,A), TU[i] for U h TA[i] for A { , TB[i] for B h TC[i] for C h and TD[i] for 



t No summation over i is implied here. 
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D{. The names of other variables used are explained in the programme. All the 
required equations are taken from the body of the paper. 

(* Input the data 

nf=number of flavours 

nc=number of colours 

d= 4 if nc is even, d= 2 if nc is odd 

r=rank=nc/2 if nc is even, r=rank=(nc-l)/2 if nc is odd 
mu= -1/2 

m[l] , m[2] , m[nf] : bare masses of the multiplets 

T[0] , T[l] , T[nf] : symmetric polynomials in the squared 

masses 

u[0], u[2], u[2r-2]: moduli 

and run *) 

(* Known moduli for SO(nc) are the highest one, which is 1, and 
all the odd ones, which vanish *) 

Do[u[j]=0, {j, 1, 2r-l, 2}] 

u[2r]=l; 

(* Definition of the curve W. For simplicity the quantum scale 
has been set to 1 *) 

p[x]=Sum[u[j] x A j , {j, 0, 2r}] ; 

W [x]=Collect [Expand [p[x]A2-xAd Product [(xA2-m[j] A2) , { j , 1, nf, 
1}]], x]; 

Print ["pM = \ p[x]]; 

Print [ M W[x] = ", W[x]] ; 

(* Discriminant of the curve *) 
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DxW=Simplif y [D [W [x] , x] ] ; 

R=Factor [Simplify [Resultant [W [x] , DxW, x]]]; 
Print ["R = " , R] 

(* Definition of the list of initial values IV [j] that will 
close the recursions *) 

Do[P[j]=IV[j] , {j, 0, 4r-2, 2}] 

initials=Table[IV[j] , { j , 0, 4r-2, 2}]; 

Print ["Initial values = ", initials] 

(* Recursion relation for the P[j] *) 

P[n_Integer]=l/(n+l-4r(2+mu)) (Sum[(-l)A(nf-i) 
(n-4r+l-(l+mu) (2i+d)) T[nf-i] P [n+d-4r+2i] , {i, 0, nf, l}] + 
2Sum[((l+mu) (j+2r)-(n-4r+l) ) u[j] P[n-2r+j], { j , 0, 2r-2, l}] + 
Sum[((l+mu) (i+j)-(n-4r+l)) u[i] u[j] P[n+i+j-4r], {i, 0, 2r-2, 1}, 
{j, 0, 2r-2, 1}]); 

(* The required values for the P[j] are computed first for the 
sake of efficiency, then expressed in terms of the initial values, 
and finally stored as the new variables PP[j] *) 

Do [PP[j]=Collect [Simplify [P [j] ] , initials], { j , 0, 8r-4, 2}] 

(* Recursion relation for the Q[j] in terms of the PP[j] *) 

Q[n_Integer]=-l/(n+l-4r(l+mu)) (Sum[(i+j-4r) u[i] u[j] PP[n+i+j], 
{i, 0, 2r-2, 1}, {j, 0, 2r-2, l}]+ Sum[(-1) A(nf-i) (4r-2i-d) 
T[nf-i] PP[n+d+2i], {i, 0, nf, l}]+ 2Sum[(j-2r) u[j] PP[n+j+2r], { j , 
0, 2r-2, 1}]); 

(* Values of P[j] or Q[j] can always be expressed as linear 
combinations of the above initial values, with coefficients that are 
polynomials in the moduli *) 



27 



(* Computation of the M matrix relating the Q[j] to the P[j] *) 

Do [QQ [ j ] =Collect [Simplify [Q [ j ] ] , initials], { j , 0, 4r-2, 2}] 

M=Table [Coefficient [QQ [i] , IV [j]], {i, 0, 4r-2, 2}, { j , 0, 4r-2, 
2}]; 

(* Next we check that the determinant of the M matrix, DM, 
contains the same factors as the discriminant R of the curve, 
possibly with different multiplicities *) 

DM=Factor [Det [M] ] ; 

Print ["DM = ", DM] 

(* Definition of modular derivatives of the Q[j] in terms of the 
P[j] *) 

DQ[n_Integer, j .Integer] =2Sum [u [i] PP[n+j+i], {i, 0, 2r, 1}]; 

(* Computation of the matrices DD[j] expressing modular 
derivatives of the Q[j] as linear combinations of the P[j] *) 

Do[ 

Do[ DQQ[n, j] =Collect [Simplify [DQ [n, j]], initials], {n, 0, 
4r-2, 2}]; 

DD[j]=Table [Coefficient [DQQ[n, j] , IV [i]], {n, 0, 4r-2, 2}, {i, 
0, 4r-2, 2}], 

{j, 0, 2r-2, 2}] 

(* Inversion of M. We pull its determinant, DM, to the left-hand 
side of the equations. SIM is the matrix of cof actors of M *) 

SIM=Simplif y [DM*Inverse [M] ] ; 

(* Computation of the matrices TU[j] of coefficients in the 
expansion of the modular derivatives of the Q[j] as linear 
combinations of the Q[j] *) 
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Do[TU[j]=Simplify[DD[j] .SIM] , { j , 0, 2r-2, 2}] 

(* Partition of the TU[j] into submatrices : TA[j] , TB[j], 
TC[j], TD[j] *) 

Do[ 

TA [ j ] =TU [j ] [ [Range [1 , r] , Range [1 , r] ] ] ; 
TB [ j ] =TU [ j ] [ [Range [1 , r] , Range [r+1 , 2r] ] ] ; 
TC [j ] =TU [j ] [ [Range [r+1 , 2r] , Range [1 , r] ] ] ; 
TD [j ] =TU [ j ] [ [Range [r+1 , 2r] , Range [r+1 , 2r] ] ] , 
{j, 0, 2r-2, 2}] 
Do[ 

Print [" Below is modulus ", u[j] , " "] ; 

Print ["TA[", j,"] = ", TA [j] ] ; 
Print [""] ; 

Print ["TB[", j,"] = ", TB [j] ] ; 
Print [""] ; 

Print ["Det[TB[", j, "]] = ", Factor [Det [TB [j] ]]] ; 
Print [""] ; 

Print ["TC[", j,"] = ", TC[j]]; 
Print [""] ; 

Print ["TD[", j,"] = ", TD [j] ] ; 

Print [""] , 

{j, 0, 2r-2, 2}] 

(* One can check that the determinants of the TB[j] contain the 
factors present in the discriminant R of the curve. *) 

29 



